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Spin rectification in a single crystal Fe/Au/Fe sandwich is electrically detected for collinear and 
non-collinear magnetization and external magnetic field configurations. The line shape, line width 
and signal polarity are analysed. The spin rectification theory has been much extended by taking 
the magneto-crystalline anisotropy and shape anisotropy into account, which explained non-collinear 
resonances and agrees very well with experimental data. Thus, a comprehensive understanding of 
spin rectihcation in ferromagnetic metal was demonstrated in this work. 
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A decade ago, spin dynamics in ferromagnetic materi¬ 
als was electrically detected via the spin diode effect in 
magneto tunnel junctions and the bolometric effect 
in thin films [3|, 1^, which triggered a rapid development. 
Later, more methods were developed, such as the spin 
pumping effect the (inverse) spin hall effect 0 
and the spin rectification effect (SRE) @,0. The SRE 
dominates the electrical voltage induced by ferromagnetic 
resonance (FMR) in a ferromagnetic metal 0 . A pre- 
cessing magnetization leads to a periodically changing 
resistance through magneto resistance. The periodically 
changing resistance couples with the microwave current 
flowing inside and generates a DC voltage, this is the 
SRE. Such a method became the most popular method 
in electrical detection of FMR because of its high sensi¬ 
tivity, simple sample structure, and experimental set up. 
It was applied to different materials and structures with 
accurate agreement between theory and experimental re¬ 
sults on both line shape and line width |1QI42Q|. Such 
line shape analysis is useful for distinguishing spin rec¬ 
tification from spin pumping and inverse spin hall effect 
00. Line width is also important for determining ad¬ 
ditional damping due to spin pumping as well as intrinsic 


Gilbert damping [5|, |10|, l2l|, |22| . All the previous studies 


of line shape and line width were performed in a collinear 
case where the magnetization is aligned parallel with the 
external magnetic field. However, in ferromagnetic thin 
film, the magnetization orientates along an effective field 
direction rather than the external magnetic field direc¬ 
tion, especially when the internal magnetic field, such 
as magnetic anisotropy field and demagnetization field, 
is comparable to the external magnetic field. In such a 
non-collinear case of the magnetization and the external 
magnetic field, the line shape and line width analyses of 
spin rectification haven’t been systematically studied yet. 

In this work, we experimentally studied the line shape 
and line width of spin rectification in a non-collinear case 



FIG. 1. RHEED patterns with the electron beam e~ along 
MgO (100) of (a) Fe (7 nm)/MgO, (b) An (4 nm)/Fe (7 nm) 
/MgO, and (c) Fe (3 nm)/Au (4 nm)/Fe (7 nm)/MgO. (d) a 
sketch of measurement geometry. 


for a sample with strong anisotropy. We also extend the 
spin rectification theory from a collinear case into a non- 
collinear case by considering all anisotropy effects. Thus, 
we present a comprehensive understanding of spin recti¬ 
fication in a metallic system. 

To achieve a system with strong anisotropy, we de¬ 
signed ultra-thin single crystal Fe/Au/Fe sandwich on 
MgO (001) substrate by molecular beam epitaxy in a ul¬ 
tra high vacuum chamber. The substrate was cleaned 
by annealing at 680°G for 45 minutes. Then, a 7-nm- 
thickness Fe was prepared at room temperature and an¬ 
nealed at 250°G for 3 minutes until the high crystalline 
quality achieved as indicated by a sharp reflection of 
high-energy electron diffraction (RHEED) pattern, as 
shown in Eig. 1 (a). A 4-nm-thickness An was then 
epitaxially deposited at room temperature. A 3-nm- 
thickness of Ee was then epitaxially deposited. Eurther, 
a 5-nm-thickness MgO layer was deposited on top for 
protection. The RHEED patterns shown in Eig. HJa)- 
(c) indicate the smoothness of each layer surface and the 
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high crystalline quality of the sample. In addition to 
the shape anisotropy, the single crystal Fe ultra-thin film 
on MgO (001) has a strong four-fold anisotropy in plane 
with the easy axis along the Fe [100] and the hard axis 
along Fe [110] [23|, and the two Fe layers with different 
thickness perform different magnetic anisotropy field 0 , 
which have all been confirmed in out measurement. Both 
magneto-crystalline anisotropy and shape anisotropy in 
Fe/Au/Fe sandwich allow us to study the non-collinear 
spin rectification in this work. 

As shown in Fig. [Il[d), the tri-layer sample was pat¬ 
terned into a strip along the Fe [100] easy axis with di¬ 
mension of 20 jumx 3 mm using standard photo lithogra¬ 
phy. A microwave was applied into the strip directly, and 
most microwave current flows inside of An layer due to 
high conductivity. Thus, the microwave magnetic field 
on the bottom layer has a phase shift of tt with that 
in the top layer. The microwave was modulated with 
a frequency of 8.33 kHz. Voltage was measured along 
the strip using lock-in technique. An external magnetic 
field H was applied to the strip with orientation defined 
in Fig. 1(d). Spin rectification voltage was measured by 
sweeping the external magnetic field at a fixed microwave 
frequency. In this work, microwave power is 100 mW. 
Before showing our experimental results, we have ruled 
out the magnetization coupling [25| and spin current cou¬ 
pling 0 between two Fe layers experimentally (experi¬ 
mental evidence was no shown in this work). Therefore, 
we treat both Fe layers independently. And, we have 
carefully checked the special condition Hi for pure spin 
pumping, and the signal is ignorable comparing to that 
of spin rectification. Thus, we were allowed to study the 
line shape, line width, and polarity of the pure spin rec¬ 
tification signal in both collinear and non-collinear cases. 

Figure [2] shows the results when H is applied near 
Fe [110] direction in the film plane, which is the hard 
axis of four-fold magneto crystalline anisotropy. Fig. 
Ela) shows a sketch of in-plane configuration measure¬ 
ment with (fH ~ 45° and 'dn = 0°. For this case, 
when H is larger than the saturation field, the magne¬ 
tization M will lie almost parallel to the H direction, 
while if H is smaller than saturation field, M will be 
pulled out of the collinear configuration, and the rela¬ 
tive angle between M and H is determined by the com¬ 
petition between the Zeeman energy and four-fold mag¬ 
neto crystalline anisotropy energy. Fig. Hb) shows uj- 
H dispersion plot, with normalized rectification voltage 
amplitude mapped into rainbow color scale as the indi¬ 
cator marks. Dispersion curves are calculated by solv¬ 
ing Landau-Lifshitz-Gilbert (LLG) equation [23|, and we 
get four-fold magnetic anisotropy field jHoHi = 0.073 T 
(black solid line) and ijlqHi = 0.026 T (grey solid line) of 
each Fe layer by fitting the measured data in Fig. [^b) . 
Due to the Fe-thickness dependence of anisotropy |24| . 
we can identify the dispersion curve traced by the black 


(^) In-plane Configuration 




(^) • 7 nm Fe FMR\\ ▼ 7 nm Fe FMR^ o 3 nm Fe FMR\\ 



FIG. 2. VsR measurement with H along hard axis Fe [110] in 
plane, (a) The sketch of the in-plane configuration measure¬ 
ment; (b) Ct;-H dispersive image plot, the grey solid line is the 
fitting curve of 3 nm Fe, and the black solid line is the fitting 
curve of 7 nm Fe; both dispersion curves have two brunches: 
brunch I is FMR\\ brunch, and brunch II is FMR]^ brunch, 
(c) Typical curves in in-plane configuration, solid circles (•) 
indicate peaks belong to FMR\\ brunch in 7 nm Fe, solid tri¬ 
angles (t) indicate peaks belong to FMR]^ in 7 nm Fe, hollow 
circles (o) indicate peaks belong to FMR\\ in 3 nm Fe. 


solid line as originating from the 7 nm Fe layer and the 
curve traced by the grey solid line as originating from 
the 3 nm Fe layer. These two dispersive curves cross 
at iiqH = ± 0.065 T, and the independence of the two 
dispersive curves near the cross indicates the magnetic 
coupling between the two FM layers is very weak. Both 
the cj -H dispersive curves have two brunches, as shown 
in Fig. [2][b). In brunch I, the resonance field increases as 
the frequency increases, here H is larger than the satu¬ 
ration field and thus M || H; we define the resonance in 
this situation as FMR\\ brunch. In brunch II, the reso¬ 
nance field decreases as the frequency increases, here H 
is smaller than the saturation field and thus M ^ H; we 
define the resonance in this situation as FMR]^ brunch. 
Fig. [2][c) shows some typical curves measured in this con¬ 
figuration at various microwave frequencies between 7.5 
GHz and 9.5 GHz. All resonance in the curves are anti¬ 
symmetric Lorenz line shape, indicating the phase shift 
between microwave field h and microwave current j is al- 
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most the integers of tt in this device [ll|. In addition to 


the rectification voltage observed at the FMR fields of the 
3 nm Fe and 7 nm Fe, a non resonant rectification signal 
is observed around /io77=0; this signal arises due to the 


shall focus our study only on the resonance rectification 
voltage. From Fig. EJc), we summarize the main fea¬ 
tures of the SRE measured in the in-plane configuration 
by the following Eqs. ([T]): (a) all voltage signals change 
their polarity when the applied magnetic field reverses; 

(b) the voltage polarity in the 7 nm Ee FMR\\ brunch is 
opposite to the polarity in the 3 nm Ee FMR\\ brunch; 

(c) the voltage polarity in FMR\\ brunch is opposite to 
the polarity in FMR]^ brunch. 


At = 0° 


V{H) = -V{-H) 

Vfct ^Fes 


I^Fe.l ■ 
Vfmr\\ 
\VfMRii 


\VFes\ 

VfmR)^ 

\VfmrA 


(la) 
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(^) Out-of-plane Configuration (b) g 



magnetic field reverses. 




|26|. In this paper we 
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Fe/Au/Fe 




PoH(T) 

(^) • 7 nm Fe FMR^\ ▼ 7 nm Fe FMR^ O 3 nm Fe FMR\\ V 7 nm Fe FMi?| 



PoH (T) 


Eq. ([Taj) is in agreement with the literatures HI, [HI , 
and Eq. ([Tb|) describes the polarity difference in two Ee 
layers due to the phase shift of the microwave magnetic 
field. Eq. ([Tc|) indicates that in the in-plane configuration 
the polarity of Vsr changes its sign for the case where 
M and H are non-collinear. And from Eig. (|2|)(c), the 
resonance peaks in FMR)^ brunch is much broader than 
in FMR\\ brunch. 


In addition to magneto anisotropy, shape anisotropy 
is also able to affect the relative angle between M and 
H. Eig. [3] shows the results when H is applied almost 
perpendicular to film plane, with Eig. [Sfa) showing a 
sketch of out-of-plane configuration measurement with 
with (fH = 0° and ~ 90°. In this configuration, 
when H is larger than the saturation field, M || H, and 
when H is smaller than the saturation field, M H. 
The relative angle between M and H is determined by 
the competition between the Zeeman energy and shape 
anisotropy energy. Eig. [3](b) shows cj-H dispersion plot, 
with normalized rectification voltage amplitude mapped 
into rainbow color scale as the indicator marks. We can 
identify the dispersion curve traced by the black solid 
line as originating from the 7 nm Ee layer and the curve 
traced by the grey solid line as originating from the 3 
nm Ee layer. Both dispersion curves also have FMR\\ 
brunch and FMR]^ brunch. Eig. He) shows several typ¬ 
ical curves measured in this configuration at various mi¬ 
crowave frequencies between 3.8 GHz and 4.2 GHz. All 
the resonance peaks show the Lorenz line shape, and we 
describe the key features by the following Eqs. (|2|): 


FIG. 3. Vsr measurement with H pointing out of the film 
plane, (a) The sketch of the out-of-plane configuration mea¬ 
surement; (b) cj-H dispersive image plot, the grey solid line 
is the fitting curve of 3 nm Fe, and the black solid line is 
the fitting curve of 7 nm Fe; both dispersion curves have two 
brunches: brunch I is FMR\\ brunch, and brunch II is FMR]^ 
brunch, (c) Typical curves in out-of-plane configuration, solid 
circles (•) indicate peaks belong to FMR\\ brunch in 7 nm Fe, 
solid triangles (t) indicate peaks belong to FMR\^ brunch in 
7 nm Fe, hollow circles (o) indicate peaks belong to FMR\\ 
brunch in 3 nm Fe, and hollow triangles (v) indicate peaks 
belong to FMR]^ brunch in 3 nm Fe. 


At ifH = 0°,79i:/ ^90° : 

V{H) = V{-H) 

Vfct _ ^Fes 
I^Fe.l " |We3| 

Vfmr^\ Vfmrj^ 

\VfmR\\\ \VfmR:^^\ 

Equations (j2|) are quite different from Eqs. ([T]). Eq. 
(|2a|) shows the voltage signal keeps the same polarity 
when H reverses, which indicates the spin pumping and 
the inverse spin hall effect is ignorable in our measure¬ 
ment Hi- Eq. (|2bj) shows the signal polarity in two Ee 
layers are the same and Eq. (|2c|) shows the signal polar¬ 
ity in FMRj^ brunch keeps the same as in FMR\\ brunch. 
Erom Eig. (|3|) (c), the resonance peaks in FMRj^ brunch 
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is also much broader than in FMR\\ brunch. Comparing 
Fig. (|2j) and (j3j), the SRE signal in FMR]^ brunch has 
the same behaviour as in FMR\\ brunch when chang¬ 
ing the measurement configuration and rf magnetic field 
direction. And comparing Vsr in two brunches, the sig¬ 
nal polarity is opposite in the in-plane configuration and 
keep the same in the out-of-plane configuration. 

So far in the literatures, SRE was systematically stud¬ 
ied only in the configuration with M || H, and the rec¬ 
tification voltage is described by a formula as a function 
of H [l^. Since M and H are non-collinear in 
brunch, the conclusions in previous studies are not suit¬ 
able here any more. However, the M alignment is always 
parallel to the effective field Hg// rather than H. Thus, 
He// instead of H should be taken into account espe¬ 
cially in ferromagnetic systems with strong anisotropy 
and demagnetization. Hg// is determined by the free 
energy F of the system. Considering the single crystal 
magnetic thin film in our case with Zeeman energy, mag¬ 
neto anisotropy energy, shape anisotropy and demagne¬ 
tization energy, one can get the free energy F and the 
effective field Heff ns follows: 

F = —jioMH [cos Oh cos Om cos((^m — ^h) + sin Oh sin Om] 

+ ^IJ.oMef/sm^0M + Ks,pCOS^9MCOS^{(fiM - <fs,p) 

+ (sin^20M + cos^0Msin^2(^M) (3a) 



d^F d^F / d'^F yi 

Kd^PMdOMj {0M vm) 
“(3b) 

Here (fM, 'Om^ and ^h are the angles of M and 
H, as defined in insets of Fig. |4l[a) and Fig. [5][a), 
/io is susceptibility in vacuum, Mgjj is effective mo¬ 
ment, Ks^p is uniaxial anisotropy constant, ips,p is the 
angle of easy axis of uniaxial anisotropy, and Ki is four¬ 
fold anisotropy constant. Putting effective field calcu¬ 
lated from Eq. m and the microwave magnetic field 
^x'Y'Z' = (O, hy/C05 ((5) O) in to LEG equation, we 
can get the dynamic magnetization m. Here b is the 
phase of the microwave field h, and in our system we 
define = 0 in 7 nm Fe and 6 = tt in 3 nm Fe. Spin rec¬ 
tification voltage is described as Vsr = {j * AR), here j 
is microwave current in the system, and AR ex Re{m) 
is resistance variation within the system due to AMR 
and spin procession. Thus we can derive the SRE in the 
in-plane configuration: 

Vsr = a * Re{xT)hY'Cos {ipM + sin (2(/9 m) (4) 

and the SRE in the out-of-plane configuration: 

Vsr = a * Re (xl) hy' sin {20m) (5) 


Mo 

{ixoMcosOMf 


with 


A 


jx' ar 
2M 


Re(xL) = 
Re( xt)= 




) 


(ooir — + 4uj‘ir 

V ^e// J J^eff 

2auj‘^ujM^^Heff 


Here jx' is the microwave current amplitude, Re(xL) 
and Re(xT) are respectively the real parts of diagonal 
and non-diagonal elements of dynamic susceptibility ten¬ 
sor, uj is the applied microwave frequency, uum = T^f, 
^Heff — l^eff^ 7 is gyromagnetic ratio, and a is damp¬ 
ing constant. As shown in Eq. (j4]) and ([5]), Vsr is a 
function of the effective field Hg// instead of the applied 
field H, thus Vsr cannot be fitted by a simple formula. 
To analysis the SRE, we first get lpm and Mm as func¬ 
tions of H by minimizing system Free energy F as shown 
in Eq. (j3aj), then calculate effective field Hg// by Eq. 
(|3b|), and finally calculate Vsr by Eq. (|4|) and ([5]). 

Figure [U shows the comparison between calculation 
and experimental results in the in-plane configuration. 
Fig. [4|(a) is a typical experimental curve measured with 
the microwave frequency of 10 GHz, and (b) shows the 
calculation curve with the microwave frequency fixed at 
10 GHz, the effective field Heff as a function of H in 
the in-plane configuration is shown in (c). Here we use 
epH = 44.6° and i^h = 0° for the in-plane configuration, 
= 0.073 T, poMeff = 1.7 T for 7 nm Fe, and 
IIqHi = 0.026 T, poMeff = 1.4 T for 3 nm Fe. These 
parameters are all determined by the dispersion curves 
in Fig. EJb). And we use a = 0.006 calculated from 
line width, and set A^^g^ = 5 x A^^gg to best repre¬ 
sent the experimental conditions. The calculation re¬ 
sults agree well with experimental results. From Eq. 
0, Vsr is determined by the real part of diagonal el¬ 
ements of dynamic susceptibility tensor Re (xl) which 
is anti Lorentz line shape, so Vsr is anti Lorentz line 
shape as shown in Fig. S^b) and confines with experi¬ 
mental result. Since Vsr oc cos {(pM + sin (2(/9 m), and 
when H reverses, the Heff and M will reverse, which 
corresponds to (pM + tt and Mm + tt, Vsr will change its 
polarity when H reverses as shown in Fig. (H^b) and con¬ 
fines with Eq. (p!a|) . And Vsr has opposite polarity in 
7 nm Fe layer and 3 nm Fe layer, as shown in Fig. |4l[b) 
and confines with Eq. (0, because in 7 nm Fe layer and 
3 nm Fe layer, the phase 6 of rf magnetic field h has a 
difference of tt. As shown in Fig. lU^c), the Heff will in¬ 
crease as H increases when M and H are collinear, which 
means spin procession is in-phase when H > Hfmr 
and out-of-phase when H < Hfmr [13 ? while the Hgj f 
will decrease as H increases when M and H are non- 
collinear, which means spin procession is in-phase when 
H < Hfmr and out-of-phase when H > Hfmr [131- 
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FIG. 4. Comparison of experiment and calculation results in 
the in-plane configuration, (a) Experiment curve with mi¬ 
crowave frequency of 10 GHz, the inset shows non-collinear 
conhguration of M and H in the in-plane conhguration. (b) 
Calculation curve with microwave frequency hxed at 10 GHz, 
the result is agree with Eqs. (Hj. Solid lines in (c) are calcu¬ 
lated effective held as a function of the applied held of 7 nm 
Ee layer and 3 nm Ee layer in the in-plane conhguration and 
the dashed line indicates the position of effective held which 
satishes resonance condition with microwave frequency of 10 
GHz. 


Near the resonance position as indicated by dashed line in 
Fig. m^c), {H - Hfmr) / {Heff — Hfmr) > 0 in FMR\\ 
brunch, while {H - Hfmr) / {Hef f - Hfmr) < 0 in 
FMRj^ brunch. And Since the sign of Vsr is deter¬ 
mined by Vsr has the opposite polarity in 

FMR\\ brunch and FMRj^ brunch when Vsr is plotted 
as a function of H, as shown in Fig. SJ^b) and conhnes 
with Eq. 

Our theory also works in the out-of-plane conhgura¬ 
tion. Fig. [5] shows the comparison between calculation 
and experimental results in the out-of-plane conhgura¬ 
tion. Fig. [Sfa) is a typical experimental curve measured 
with the microwave frequency of 4 GHz, and (b) shows 
the calculation curve with the microwave frequency hxed 
at 4 GHz, the effective held He// as a function of H in 
the out-of-plane conhguration is shown in (c). In calcula¬ 
tion, we use (fH = 0° and = 89.4° for the out-of-plane 


conhguration, and keep the other parameters the same as 
those used in the in-plane conhguration. From Eq. m, 
Vsr is determined by the real part of non-diagonal ele¬ 
ments '^f dA/^namir* fpn«nr A/n-.^ AArhiph fg 

• 7 nm Fe FMR\\ ▼ 7 nm Fe O 3 nm Fe FMR\\ V 7 nm Fe 



- 2-10 1 2 
goH (T) 


EIG. 5. Gomparison of experiment and calculation results 
in the out-of-plane conhguration. (a) Experiment curve with 
microwave frequency of 4 GHz, the inset shows non-collinear 
conhguration of M and H in the out-of-plane conhguration. 
(b) Galculation curve with microwave frequency hxed at 4 
GHz, the result is agree with Eqs. Solid lines in (c) are 
calculated effective held as a function of the applied held of 7 
nm Ee layer and 3 nm Ee layer in the out-of-plane conhgura¬ 
tion and the dashed line indicates the position of effective held 
which satishes resonance condition with microwave frequency 
of 4 GHz. 

Lorentz line shape, as shown in Fig. ©(b), and conhnes 
with experimental results. Since Vsr oc sin( 26 >M), Vsr 
keeps the same polarity when H reverses (conhnes with 
Eq. (j2aj)), and keeps the same polarity in 7 nm Fe and 
3 nm Fe layer (conhnes with Eq. And Since the 

sign of Vsr is determined by ujReff ^ ^SR polarity keeps 
the same in FMR\\ brunch and FMR]^ brunch. 

The calculation and the experimental results of the 
SRE in 7 nm Fe layer are listed in Table HI Our theory 
well describes the line shape and polarity of the SRE in 
the general conhguration with M and H. Also our theory 
conhrms the broaden of linewidth AiT" when M and H 
are non-collinear qualitatively. However, the broaden of 
linewidth in experiment is larger, and the quantitative 
analysis still needs further discussions. 
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Measurement Configuration 

Line shape 

polarity 

rqHemr (T) rqAH (T) 


(fH 

= 44.6°, = 0°, / = 10 GHz 

Anti-Lorentz( Exp) 

— (Exp) 

0.14(£;a,p) 

0.0053(£,p) 

M II H 

Anii-LoTeniz (^Cai) 

— (Cal) 

0.14(c'aZ) 

0.0024(ca,) 


= 0°,^?/f = 89.4°,/ = 4GHz 

Loientzi^Exp) 

~^(Exp) 

1.96(£;a,p) 

0.0095(£,p) 


Lorentz(c'aZ) 

+ (CaZ) 

l-97(CaZ) 

0.0072(cai) 


(fH 

= 44.6°,= 0°,/ = 10 GHz 

Anti-Lorentz(£;a,p) 

-^(Exp) 

0.027(^Exp) 

0.010(£;2;p) 

M It H 

Anii-LoTeniz (^Cai) 

+ (CaZ) 

0.026(c'aZ) 

0.0036(coi) 


= 0°,i?/f =89.4°,/ = 4GHz 

Lorentzi^Exp) 

-^(Exp) 

^ •l^(Exp) 

0.054^£;3.p, 


Lorentz(c'aZ) 

+ (CaZ) 

l-80(CaZ) 

0.013(Cal) 


TABLE L The calculation and the experimental results of the SRE in the M || H and the M H configuration in 7 nm Ee 
layer with different measurement geometry. The positive polarity of the SRE is defined as Vsr/\Vsr\ > 0 when H < Hrmr- 
The subscript Exp indicates the result is extracted from the experimental data, and the subscript Cal indicates the result is 
extracted from the calculation data. 


In conclusion, we studied Spin Rectification Effect in 
an epitaxial Fe/Au/Fe tri-layer system with strong mag¬ 
neto anisotropy and shape anisotropy. In addition to the 
SRE when M and H are collinear, we study Vsr for the 
case where M and H are non-collinear. The different be¬ 
haviour of Vsr in different configuration of M and H are 
due to the different relationship of the Hg// depending 
on H. By considering Hg// instead of H in ferromag¬ 
netic system, we extend the SRE theory for all M and 
H configurations in different measurement configuration. 
These equations will help further understanding of spin 
transport in ferromagnetic systems, especially when M 
is not parallel to H. 
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